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We present the non-relativistic QCD (NRQCD) prediction for the polarization of the J/ip pro- 
duced in 6 ^ J/tp + X, as weU as the hehcity-summed production rate. We propose that these 
observables provide a means of measuring the three most important color-octet NRQCD matrix 
elements involved in J/ip production. Anticipating the measurement of the polarization parameter 
a, we determine its expected range given current experimental bounds on the color-octet matrix 
elements. 



I. INTRODUCTION 

A rigorous theoretical framework within which quarkonium production can be studied is provided by the nonrel- 
ativistic QCD (NRQCD) factorization formalism developed by Bodwin, Braaten, and Lepage This approach is 
based on NRQCD an effective field theory that can be made equivalent to full QCD to any desired order in 
heavy-quark relative velocity. The NRQCD factorization formalism is not a model, but rather a rigorous derivation 
within NRQCD of a factorized form for quarkonium production and decay rates. 

A central result of the NRQCD factorization formalism is that inclusive quarkonium production cross sections 
must have the form of a sum of products of short-distance coefRcients and NRQCD matrix elements. The short- 
distance coefficients are associated with the production of a heavy quark-antiquark pair in a specific color and angular- 
momentum state. They can be calculated using ordinary perturbative techniques. As to the NRQCD matrix elements, 
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they parameterize the effect of long-distance physics such as the hadronization of the quark- antiquark pair. These 
can be determined phenomenologically. 

The power of the NRQCD formahsm stems from the fact that factorization formulas for observables are essentially 
expansions in the small parameter , where v is the average relative velocity of the heavy quark and anti-quark in the 
quarkonium bound state, ~ 0.3 for charmonium, and 0.1 for bottomonium. NRQCD w-scaling rules |3j allow us to 
estimate the relative sizes of various NRQCD matrix elements. This information, along with the dependence of the 
short-distance coefficients on coupling constants, permits us to decide which terms must be retained in expressions for 
observables so as to reach a given level of accuracy. Generally, to leading order, factorization formulas involve only a 
few matrix elements, so several observables can be related by a small number of parameters. 

NRQCD predictions therefore presuppose QCD, factorization, and a reasonably convergent v'^ expansion. If the 
quantitative predictions of NRQCD fail, one of these three assumptions is failing, most likely the validity of the 
expansion. 

In many instances of J/tp production, the most important NRQCD matrix elements are (Of (■^S'l)), {Ofi^Si)), 
{0^(^80)), and {Of{^Pj)). These four quantities parametrize the hadronization into a J/ip boundstate of a cc 
pair produced initially with the stated quantum numbers (angular momentum ^^'^^Lj, color quantum number 1 or 
8) . Previous to the development of the factorization formalism of Ref . , most J/ ip production calculations took 
into account only the hadronization of cc pairs initially produced in the color-singlet state, as parameterized by 

Recently, in regard to the three most important color-octet matrix elements involved in J/ip production, it has 
been noticed that there prevails a certain inconsistency between the values measured at CDF fl and those measured 
in other processes . We propose here that inclusive J/ 4' production from 6-decay can be of avail in unravelling this 
difhculty by serving as a supplementary arena for measuring the contentious matrix elements. 

Inclusive production of J/ip from fe-decay provides two measurable combinations of the matrix elements. The first 
one is the helicity-summed rate r(6 — > J/ip + X). The second combination concerns the polarization parameter a 
appearing in the electromagnetic decay rate of J/ip to lepton pairs: 

-U ^ n+fi-{0)) ocl + acos'^e , (!) 

acost) 

where the polar angle 6 is defined in the J/^ rest frame for which the z-axis is aligned with the direction of motion 
of the J/ip in the lab. 

The short-distance physics for the process b — > J/ip + ^ is described by the four-quark Fermi interactions b ccs 
and b —^ ccd, at the mb scale. The Feynman diagram is shown in Fig. 1 . Let us define P to be the total four- momentum 
of the cc system. The s and d quarks are assumed massless. 



II. FACTORIZATION FORMULA FOR THE SPIN-SUMMED PRODUCTION RATE 



We first calculate the factorization formula for the spin-summed J/t/i production rate. Let us define q as the relative 
momentum of the c and c quarks evaluated in the rest frame of the cc pair. The NRQCD matching formalism requires 
only that we know the algebraic form of the amplitude for small values of q. This is because the hadronization of the 
heavy quark pair into quarkonium occurs with nonnegligible probability only when |q| rric. Moreover, this is the 
regime in which the NRQCD effective lagrangian is valid. Therefore, we expand the Feynman amplitude in powers of 
q and keep only pieces up to linear order. The amplitude for b ccs is given by 



A^(cr,r;c,d,e,/;q) 
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where ^ and 77 are non-relativistic two-spinors for the c and c quarks respectively; c, d, e, / are color indices and a and 
T are the individual spins of the charm quarks. The factors involving C+ and C_ are Wilson coefficients which govern 
the scale-evolution of the four-quark Fermi interaction. (defined in detail in Ref. |^) is the Lorentz boost matrix 
that takes a three-vector from the cc rest frame to the lab frame . In the above equation, the term with the color 
Kronecker deltas contributes to the production of cc pairs in color-singlet states, and the term with color matrices 
contributes to the production of cc pairs in color-octet states. 

From Eq. ^, using the matching procedure presented in Ref [Q, we determine the spin-summed rate for b J/tp + X 
in the NRQCD factorization formalism: 
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where, upon summing over the hght quarks s and d, we have used IKbl^csP + l^cbV^'^P ~ |KbP- We have taken the 
inital 6-quark to be unpolarized. Here (Ojf) = (OjO^lO) are NRQCD J/tp production matrix elements. We consider 
only the leading color-singlet piece and the leading color-octet pieces in the relativistic w^-expansion. Our calculation 
of the coefficients in the above spin-summed factorization formula concurs with the results presented in Ref . 1^ . 

According to the NRQCD ti-scaling rules, the color-octet matrix elements in Eq. |^ are all expected to be suppressed 
by v'^ with respect to the basic color-singlet matrix element {Of {^Si)). Thus, at the outset, one would not expect 
the octet contributions to play a major role in the production of J/t/j articles. However, as pointed out in Ref. 
the color-singlet coefficient 2C+ — C_ decreases significantly as it is evolved down from its value of 1 at the scale Mw 
to its value of roughly 0.40 at the scale rrib. On the other hand, the color-octet coefficient (C+ + C-)/2 increases 
slightly from 1 at the scale AIw to roughly 1.10 at the scale mf,. For this reason, the octet contributions are actually 
as important as, or more important than, the basic singlet contribution. In fact, the short-distance coefficients in the 
color-octet terms are some 50 times larger than those in the singlet term! 



III. FACTORIZATION FORMULA FOR HELICITY RATES 



In the previous section, we discussed how the serendipitous Wilson-coefficient enhancement of the octet contributions 
enables us to measure the color-octet matrix elements contributing to the helicity-summed rate r(6 — > J/ip + X). This 
enhancement can be further exploited to determine the color-octet matrix elements by considering the polarization 
of the produced J/ip's, as measured by the parameter a in Eq. |l| . It can be calculated via the formula 

a(+)+aH-MO) 

ai+) + ai-) + 2a{0) ^' 

where a{X) represents the production rate of J/tp with helicity A. 

Since a is a ratio, we need only calculate the relative sizes of the various contributions to the helicity rates. Thus, we 
need only calculate the effective Feynman amplitude squared, evaluated in the b rest frame, with the J/ip momentum 
P oriented in the positive z-direction. In this case, the angular momentum of the J/ip in the canonical z-direction 
corresponds simply to the particle's helicity. 

Beneke and Rothstein pO) and Braaten and Chen |^ have developed techniques for deriving the production rates 
of quarkonia with specified helicities. Applying their methods to the amplitude given in Eq. ^, we obtain 

r{b ^ J/tPiX) -f X) cx 2(2C+ - C-)^{OtCSi)) [{ml - 4ml)Sxo + ^ml{l - A)] 
+ 3(C+ + C-)^{OtCS,)) [{ml - Aml)6xo + 4m2(l - A)] 

+ 9(C+ + C_)2i^»i^^ [{ml - 4m2)(l - 5xo) + 8ml - AmlX] . (5) 

In displaying Eq. |^, we have adopted the standard practice of rewriting the matrix element (08(''^'i)) using the 
relation 

{OtePj))^i2J+l){OtCPo)). (6) 

This latter relation is due to heavy-quark spin symmetry, and is valid up to relative order w^. It must be kept in mind, 
however, that 6-decay does not actually produce cc states in the ^Pq configuration at leading order in the coupling 
constants. 

Some consequences of angular momentum conservation and of the Icft-handedness of the charged-current Fermi 
interaction are illustrated in Fig. 2. The thin arrow labelled P represents the trajectory of the center of mass of 
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the cc system. The other thin arrow represents the trajectory of the s or d quark. Thick arrows denote intrinsic 
angular momenta, with Ace denoting the heUcity of the cc system prior to hadronization. Since the charged-current 
Fermi interaction couples only the left-handed parts of the fermion fields, the masslessness of the s and d quarks 
ensures that they are emitted with negative helicity. Therefore, due to angular momentum conservation, the helicity 
outcome Xcs = +1 is not allowed. This has certain obvious consequences for the J/ip helicity production rates: 
the basic color-singlet contribution (parametrized by {Of{^Si)) and involving the direct hadronization of a color- 
singlet ^Si state into a J/tp) is zero for A = -1-1, at lowest order in the relativistic expansion. The same is also 
true of the (Og ('^S'l)) contribution, which involves the hadronization of a color-octet ■^S'l state into a J/tp via the 
heavy-quark-spin-preserving L = emission or absorption of two soft gluons . 

The color-singlet contribution to the helicity rates of the J/ tp produced in 6-decay was calculated in the color-singlet 
model by M. Wise in Ref. 1 12| 16 years ago. This author presented expressions for the production rates of longitudinal 
and transverse helicities, given in terms of the color-singlet radial wave function at the origin. The result in Ref. |l^ ] 
concurs with the first line of our Eq. |5[ 



IV. PREDICTION FOR THE POLARIZATION PARAMETER ALPHA 

While experimental determinations of helicity-summed T{b J/ip + X) have already been carried out |Q, a 
measurement of the polarization parameter a is not yet available. Anticipating the availability of this latter measure- 
ment, it is interesting to determine the range of a which is predicted by the NRQCD factorization formalism. Using 
Eqs. I and |, we can express a in terms of the NRQCD matrix elements. Before preceding, however, we must first 
decide which value of m?, to use. Since we are presenting a leading-order calculation, we may in principle choose to 
use either the pole mass or the MS mass for mb, the difference between these two being merely a higher-order-in-aj 
effect. However, it turns out that the leading-order NRQCD prediction of the polarization parameter a depends quite 
strongly on the choice of mf,. Therefore we report results for a wide range of rrif,, from 4.1 GeV to 5.3 GeV. This 
range includes the values of rrib — 4.3 ± 0.2 GeV corresponding to the MS mass, and the values rrn, ~ 5.0 ± 0.2 GeV 
corresponding to the pole mass determined on the lattice p4[ . Our range is centered around mt, = 4.7 GeV, which is 
the central value for the pole mass |]l3[. This choice is consistent with the fact that the numerical value that we use 
for {Ofi^Si)) is taken fr om Ref. |14|, in which the pole mass was used in the extraction of the matrix element value. 

For concreteness, we present in the text below only our results for the choice m?, = 4.7 GeV. Using C+(TOf,) = 0.868, 
C-{mb) = 1.329 and iric = 1-55 GeV, one obtains the leading-order NRQCD prediction for the polarization parameter 

-0.39(of - 17{0t{'S,)) + 52{ot{'Po))/ml 



(Of (3^1)) + U{Ot{^Si)) + 6l{Ot{^So)) + 21l(0|(3Po))/™2 

It must be noted that a depends only weakly on {Of{^Si)) and most strongly on {Og Pq)) /ml. 

We now precede to determine the range of a which is consistent with existing information on the matrix elements 
{Of{^Si)), (Og (35*1)), (Og (^S'o)), and {Of{^Po))/Tnl and with various constraints on linear combinations of these 
quantities. We first review this information: 

• Bodwin et al. Q have determined the color-singlet matrix element (Of (^S'l)) from decay rates of J/tp to lepton 
pairs to be (Ofi^Si)) = 1.1 ± 0.1 GeV^ The error reflects theoretical uncertainties due to higher order as and 
v'^ corrections. This phenomenologically extracted value for the color-singlet matrix element is in agreement 
with the lattice calculation result also presented in Ref. p^ . 

• A constraint on the octet matrix elements is given by the requirement that the theoretical prediction for the 
spin-summed production rate T{b —> J/tp + X) given in Eq. be consistent with the experimentally measured 
value of 3.42 /iteV. This constraint can be expressed as 

0.342 GeV^ = 0.096(Of f^i)) -t- A.2l{otCSi)) + 6.76{OtCSo)) + 25.3 . (8) 

To compute the above experimental value of r(6 J/ip + X) = 3.42/ieV, we have used 

i3i?(i3(charge not determined) — > J/'0(direct) -I- X) =(0.80±0.08)%, and T(averaged over B hadrons)=1.54 pi- 
coseconds. To compute the short-distance coefficients on the right-hand-side of Eq. ||, we have taken 
Gf = 1.1664 X 10"^ GeV"2 Q ^nd \Vcb\ = 0.0381 ± 0.0021 One can attribute an error of roughly 

30% to the above equation, to reflect relativistic corrections and uncertainties in the scale used to evaluate the 
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Wilson coefficients [^|6j. It must be pointed out that the color-singlet model prediction for spin-summed J/4> 
production — obtained by setting all color-octet matrix elements to zero, and using (Of(35i)) = l.liO.lGeV^ 
— is roughly one-third of the experimental value. However the NRQCD prediction for the &-decay width can 
be made to agree with the experimental measurements for reasonable values of the color-octet matrix elements. 

• Cho and Leibovich Q have performed a fit to CDF data over low and high pT ranges and have found 

{oti^Si)) = 0.0066 ± 0.0021 GeV^ (9) 



and 



{OtCSo)) + si^iil^^ ^ 0.066 ± 0.015 GeV^ . (10) 



• From photoproduction data, Amundson ct al. |17|] have determined 

m 



(OfCSo)) + 7 ^ ^\ ^" = 0.020 ± 0.001 GeV^ . (11) 



Note that the errors quoted in Eqs. ^, |l^, and 111] are statistical only, and that the analyses in |^ and |T^] 
did not take into account theoretical uncertainties due to next-to- leading order (NLO) corrections. However 
we expect NLO corrections to be very large (see for a discussion of NLO corrections in hadronic J/ip 
production calculations). Moreover the photoproduction results of Ref. [0 may also have significant higher- 
twist corrections pO[ . 

• Finally, due to w-scaling arguments, we expect the color-octet matrix elements to be roughly of order 
v^Oti^S,)). 

With this information in mind, we now determine the range of values of a that we can expect. Our method of 
determining this range consists of "scanning" through the three-dimensional parameter space of the color-octet matrix 
elements, and determining the maximum and minimum value of a that occur in the allowed volume. Our scan is 
subject to the following constraints: we allow {O^i^Si}) to vary in the range 

{Oti^Si)) e [0.003, .014] GeV^ ; (12) 

we allow the photoproduction-related linear combination of {0^(}Sq)) and {0^{^Pq)), to vary in the range 

(oti^S^)) + 7^^!^ e [0.0, 0.04] GeV^ , (13) 

where the range we have chosen takes into account theoretical uncertainties due to NLO corrections and possibly 
significant hight-twist corrections; we allow the b-decay-related linear combination given in Eq. ||to vary ±30% about 
the central value of 0.342 GeV^, i.e. 

0.096(of (35i)) + 4.2l(0^(35i)) + &.7&{otCSa)) + 25.3^^^^^^ e [0.24, 0.45] GeV^ , (14) 

where the range we have chosen takes into account theoretical uncertainties due to relativistic corrections as 
mentioned previously. Note that the above constraints are sufficient to insure that the helicity rates in Eq. ^ are 
positive for each value of the helicity A. 

To the above experimental constraints on the allowed parameter space, we must of course add the theoretical 
requirement that the absolute values of the octet matrix elements respect w-scaling rules, i.e. that they must not 
be unreasonably larger than u'*(Of (^5i)) ~ 0.1. Moreover, we must impose on our allowed volume of parameter 
space that the matrix element (Og (^6*0)) always be positive. No similar constraint need to be applied to the matrix 
element (Of (^Pq)). Indeed, it can be negative, for the following reason. The bare matrix element (which is of course 
positive definite) contains a power ultraviolet divergence. Such a power ultraviolet divergence must be proportional 
to a matrix element of an operator of lower dimension. In the case of {Of {'^ Pq)) /ml, the divergence is proportional 
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to (Ofi^Si)) §. This diver gence must be subtracted to obtain the renormahzed matrix element. Since the piece 
subtracted is comparable in magnitude to the bare matrix element, the difference can be negative. On the other 
hand, since the operator Of{^So) belongs to the set of lowest dimension NRQCD four-fermion operators, the bare 
matrix element {Of{^So)) cannot have power divergences. Thus the subtractions involved in renormalization cannot 
transform the positive definite bare matrix element into a negative quantity p9| ]. Due to these considerations, we 
therefore impose the additional theoretical constraint 

{OtCSo))>0. (15) 

The constraints expressed in eqs. |l^, |lj and |l5| together define the volume in the three-dimensional parameter 
space that is allowed by current experimental information and theoretical considerations. 

In defining the allowed parameter space, we have not used an experimental constraint on the CDF-related linear 
combination presented in Eq. ^ Interestingly, we find that with the imposition of all the above constraints (Eqs. |l^ 
through , this combination is limited to take values in the range 

(0^(^^o))+3 ^ °" e [0.01,0.06] GeV^ . (16) 

The central value of 0.066 quoted in En. ^ falls just outside the upper end of this range, indicating a possible 
inconsistency between the result of Ref. Q and the results from other processes considered in the present analysis. 

We now present our main result, which is the expected range of a: The maximum value for a (—0.09) is ob- 
tained when {O'^i^Si)) is at the minimum of its allowed range and when the photoproduction- related combination 
(Og (^S'o)) -l- 7 {O^ Pq)) / rn^ is at the maximum of its allowed range. The minimum value for a (—0.33) occurs in the 
opposite situation. 

So far in our discussion, we have considered only the particular choice mf, — 4.7 GeV. Table I shows how our 
findings depend on mf,. The main results are the NRQCD predictions for the range of a. Alongside the full NRQCD 
results, we display the color-singlet model predictions for a. The error bars on the color-singlet-model predictions 
reflect v'^ relativistic corrections; the relativistic corrections to the color-singlet result are of order v*, not as one 
might expect a priori, because the v'^ corrections to the helicity rates factor and cancel in the ratio a. In the other 
three columns of Table I, we treat the quantities appearing at the top of the columns as functions of the variables 
through which we scan; we report the ranges of values of these quantities that occur for the allowed parameter space. 

Our leading-order results depend strongly on the values taken for rrib and rric- As stated previously we have used 
rric = 1.55 GeV to generate the results given in Table I. For some optional choices of TOc we obtain the following 
ranges of a. Choosing mf, — 4.7 GeV and nic — 1.25 GeV we find that a E [—.53, —.12]. Choosing mi, — 4.7 GeV and 
TOc = 1.85 GeV we find that a £ [-.22, -.13]. 

In general for higher values of rub the predicted range for a enlarges. For increasing values of rric the minimum 
value we predict for a becomes greater, while the maximum value remains roughly unchanged. 

A general conclusion that can be drawn by looking at Table I is that the inclusion of octet matrix elements raises 
the predicted range for a significantly from the color-singlet prediction. 



V. CONCLUSION 



We have proposed that a measurement of the polarization of J/^p particles produced in the process b J/ip + X 
can furnish a useful new means of determining the color-octet matrix elements involved in J/^ production. This 
measurement, which is not yet available, will supplement existing experimental information, which includes extractions 
from CDF, photoproduction, and the helicity-summed rate r(6 J/ijj + X). Further experimental information would 
be of great utility, since there currently prevails a certain inconsistency between the values of the matrix elements 
measured at CDF and those values extracted from other processes. 

We have presented the factorization formula for the helicity-summed rate in Eq. ||. Rates for specified J/ tp helicity 
have been calculated using the methods of Beneke and Rothstein [|l^ and Braatcn and Chen [^ , and are presented 
in Eq. ^. The resulting expression for the polarization parameter a is given in Eq. |^ (for the choice nib = 4.7 GeV). 

In anticipation of the measurement of the polarization parameter a, we have found it interesting to determine 
the range of a which is predicted by the NRQCD factorization formalism, given current experimental information 
on the NRQCD matrix elements. We find that, for nib = 4.7 GeV, a is expected to range from —0.33 to —0.09. 
Our method of determining this range consists of "scanning" through the three-dimensional parameter space of the 
color-octet matrix elements that is allowed by the constraints expressed in Eqs. O, lO, 11^ and ll^. Upon inspecting the 
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rrib 
(GeV) 


color-singlet 

model predictions 
for a 


NRQCD 

factorization formalism 
predictions for a 


range of 
(GeV^) 


range of 
(GeV^) 


range of 

(O8*('-P0)) 

m- 

(GeV^) 


4.1 


-.27± .03 


[-.23, -.14] 


[.08,.2] 


[.1,.4] 


[-.06, -.009] 


4.4 


-.34± .03 


[-.28, -.11] 


[.03,.l] 


[.02,.2] 


[-.03, .003] 


4.7 


-.40± .04 


[-.35, -.08] 


[.009, .06] 


[0,.l] 


[-.01, .006] 


5.0 


-.45± .05 


[-.42, -.09] 


[.002, .03] 


[0,.06] 


[-.008, .006] 


5.3 


-.49± .05 


[-.47, -.12] 


[.001, .02] 


[0,.03] 


[-.005, .005] 



TABLE I. Theoretical leading-order predictions for the expected range of a. We present predictions for a in the color-singlet 
model and in the NRQCD factorization formalism. The ranges are determined by "scanning" through the allowed volume in 
the three-dimensional parameter space constrained by Eqs. (12), (13), (14) and (15) and by finding within that volume the 
maximum and minimum occuring values of a. Results are given for various values of m\,. Also given, in the last three columns, 
are the ranges of the CDF-related combination and of the two matrix elements indicated, resulting from the set of constraints. 



parameter-space volume so defined, we find that the linear combination of matrix elements that has been measured 
at CDF [|, namely (©^(^S'o)) + ?,{Ot{^ Pq)) /ml is limited to take values from 0.01 to 0.06 GeV^, a situation which 
is at the threshold of incompatibility with the value of 0.066 ± 0.015 reported in Ref. [Q. 

In Table 1, we have shown how our leading order predictions depend on mf,. In general, it can be concluded that 
the inclusion of octet matrix elements raises the predicted range for a significantly from the color-singlet prediction, 
and that the range of a broadens as mf, is increased. 

Unfortunately, due to the large number of poorly determined parameters, including rrib and mc, the expected range 
of a as predicted by NRQCD is large. 

It may be that an accurate NRQCD prediction of the polarization parameter (beyond leading order) will not 
be possible for a long time. This possibility, however, does not degrade the importance of our calculation, since an 
experimental determination of a — in conjuction with our results — will most certainly serve to tighten the constraints 
on the possible values of the color-octet matrix elements. In fact, the NRQCD prediction for a is very sensitive to 
the values of the color-octet matrix elements, and this offers the hope that a measurement of the polarization of J /tp 
produced in h -> J/V' -l- X will be instrumental in determining (©^'(^Si)), {0^g{^So)), and {Of{^Po)). 
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Figure Captions 



Figure 1. Feynman diagram for the four-quark Fermi interactions b ccs and b — > ccd. This interaction describes 
the short-distance physics for the process b J/4' + X. P is the total four- momentum of the cc system, and s, d 
and b arc the four-momenta of the corresponding quarks. The cc pair is produced initially with quantum numbers 
a, ^^~^^Lj with a = 1 or 8. It then hadronizes into a J/tp particle with helicity A. 

Figure 2. Angular momentum conservation in 6 ^ ccs and b ccd. Because the charged-current Fermi-interaction 
couples only the left-handed parts of the fermion fields, the s and d quarks (being massless) are emitted with negative 
helicity. Due to angular momentum conservation, therefore, the outcome Ace = -1-1 is not allowed. As a result, the 
(of (^S'l)) and {Og{^Si)) contributions to the J/tjj helicity rate are zero for A = at lowest order in v"^. 
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